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EARLY TIME COUPLING STUDIES USING A ID HYBRID CODE 


I. INTRODUCTION 

The understanding of atmospheric effects generated by a HANE is of 
vital importance to DNA. It is well known that HANE's cause severe 

disturbances in the atmosphere and ionosphere which impact C^I systems. In 
particular, long lasting, large scale ionization irregularities degrade 
radar and communications signals. To develop predictive capability for 
such HANE effects, DNA sponsored a major research effort at NRL leading to 
theoretical and computational models of HANE. Significant progress was 

made in this area during the past decade leading to the development of 
sophisticated raultifluid codes which incorporate in a self consistent 
manner anomalous transport due to instabilities. Recent advances in both 
plasma theory and numerical techniques allow us today to produce 
substantially more refined and accurate models of the early time HANE 
phenomenology. These advances indicate the need for a better incorporation 
of kinetic effects due to strong deviations of the particle distributions 
from Maxwellian. The new codes allow us to compute the laminar and 

turbulent dynamic ion response over long time and space scales and for a 

mixture of ionized species. In this way laminar, electrostatic, Larmor and 
turbulent coupling processes can be followed dynamically and their relative 
strength and dynamic interplay can be assessed. This paper is the first to 
describe this series of investigations, focuses on a particular process, 
early time debris-air coupling, and stresses quantitative understanding of 
the underlying physics. 

Early time coupling had been a controversial subject in the DNA HANE 
community in the I970's. Two "opposing" views were maintained: NRL, LASL 
and ARA advocated short-range coupling associated with plasma turbulence, 
while MRC argued for "Larmor coupling" associated with the ion gyration in 
the ambient magnetic field. To a certain extent, the issue was not whether 
or not these coupling mechanisms worked, but which mechanism was 
dominant. The turbulent coupling process developed at NRL was studied by 
incorporating theoretical estimates of anomalous transport coefficients in 
a ID multi-fluid code. Although this procedure yielded results consistent 
Manuscript approved February 14,1985. 








with HANE data and provided strong support for turbulent coupling, it was 
not a self-consistent technique. Moreover, it was not capable of handling 
kinetic effects such as caused by reflected Ions; such reflected ions are 
known to occur in high Mach number shocks. The present study, using a ID 
hybrid code, overcomes these deficiencies. Namely, it allows instabilities 
to be excited which affect the plasma evolution in a self-consistent 
manner, i.e., no model anomalous transport coefficients are used, and 
reflected ions are permitted. Furthermore, the ID hybrid code inherently 
includes "Larraor coupling" so that a direct comparison of the relevant 
importance and interplay of the coupling mechanisms can be made. 

The organization of the paper is as follows. In the next section we 
present an overview of the basic physical processes occurring as debris 
streams through the air. In Section III we present results of two 
simulations which demonstrate the various aspects of coupling. Finally, in 
the last section we summarize our results and discuss their implications to 
HANE phenomenology. 

II. OVERVIEW OF PHYSICAL PROCESSES 

The evolution of the debris-air plasma is characterized by three 
temporal stages: (1) magnetic field compression; (2) piston formation; and 
(3) shock formation and evolution. The bulk of this paper describes in 
detail piston dynamics; for completeness we briefly discuss the other two 
stages. 

In the initial stage, the debris streams outward from the burst point 
and picks up air electrons which merge with the debris electrons to form a 
single electron distribution. The compression of the air electrons leads 
to a compression of the magnetic field since the air electrons are "tied" 
to the field lines due to the Low value of resistivity. This is the first 
stage of the debris-air interaction and has been well-documented 
theoretically (Longraire, 1963; Sloan, 1970; Lampe and Hernandez, 1972), 
numerically (Clark et al., 1974), and experimentally (Ripin et al., 1984). 

The next stage involves the dynamics of piston formation and will be 
discussed in detail in the following section. Briefly, this stage is 
dominated by the interaction of the two ion streams, i.e., the debris ions 
and the air ions. It is this phase where various coupling mechanisms 





become apparent. The debris ions slow down, while the air ions are picked 
up; there is momentum exchange between the ion species. There is also a 
conversion of the directed kinetic energy of the debris to thermal energy 
of both ions and electrons, i.e., heating and acceleration occurs. And 
finally, an electrostatic potential forms at the leading edge of the 
magnetic field compression which causes some debris ions to be accelerated 
ahead of the magnetic piston and some air ions to be reflected. This final 
point, although well-known, can be very important to the formation of 

debris patches in the conjugate regions since it results in very high 
velocity ions. 

Lastly, the final stage involves the formation of a shock wave which 

"runs" ahead of the piston. We do not describe this phase in detail in 

this paper, but do show that a minimum time is required for high Mach 

number shock formation. This minimum time is basically the time needed for 
the shock to "run" ahead of the accelerated debris or reflected air ions. 


III. SIMULATIONS AND RESULTS 

Prior to discussing the details of the simulations, we present a 
simple model equation that elucidates the various coupling mechanisms to be 
studied. Moreover, it highlights the distinction between previous models 
(e.g., KLYSMA) and the present work. We consider the following momentum 
equation in the x direction for an air ion (i.e., radial direction 
perpendicular to the ambient magnetic field B_* Be ). 


dV 

ax 

dt 


eE„ 


*(v 


( 1 ) 


where the subscript a refers to air, fi a = Z a eB/m a c, v* is an anomalous ion- 
ion collision frequency, = ^d e x * s t * le streara:1 - n 8 debris velocity, e is 
the charge, m is the mass, and Z is the charge state. The first term on 
the RHS of (l) arises from a laminar electric field usually associated with 
the leading edge of the magnetic compression. It acts to accelerate debris 
ions and to reflect air ions. The second term in (1) is the magnetic force 
which is associated with Larraor coupling. The final term in (1) 
corresponds to turbulent "pick up" of the air and arises because of plasma 
instabilities. For plasma turbulence such that v* > It is clear that 
turbulent coupling can dominate over Larmor coupling. 








Notice that the force ft V associated 
a ay 

proportional to the value of V. 


is V a . 


ay 


0. 


with Larmor coupling is 
a y which, before any form of coupling occurs, 
Therefore the laminar force eE /m and the turbulent 


force v*|V ax - V^l ~ v*V^ will dominate initially even for v* « ft a (i. 
as long as v* > ^ a (V a y/V^)). The time evolution of V a ^ will be given by 

dV eE eV.B 

= —1 . — L. = ft V 
dt m me ad 


( 2 ) 


In deriving eq. (la) we assumed an infinitely conducting (n 0), low 8 
(8 << l) situation so that E y = V^B/c. From eqs. (1) and (2 ) it is 
obvious that the time scale for significant Larmor coupling is ft & t > 1, 
i.e., the value of V a ^ should approach Vj. At earlier times the laminar 
and turbulent forces will dominate. The turbulent forces will also 
dominate even at t » ft” 1 as long v* > ft a |V ax - Another important 

point follows from (1) and (2) if we neglect the laminar and turbulent 
forces. If L is the si 2 e of magnetic field compression then as long as V ax 
< Vj the maximum value of V a y will be V g y = ft a L. As mentioned previously 
the Larmor momentum coupling force in the x direction becomes significant 


an extremely broad magnetic field profile will be required. We will 

further comment on this point later on. 

In previous work, theoretical estimates of v* were incorporated in a 
multi-fluid code (KLYSMA) which solved momentum equations similar to (1) 
and (2) (Clark et al., 1974). However, the turbulent interaction involves 
ions interacting with fluctuating electric fields, i.e., v* ■ v*(5E). The 
present simulations incorporate such effects self-consistently. The air 
momentum equation solved, analogous to (1), is 


dV 


0B 


6E 


: 1 + ft v 

x' a ay 


(3) 


where SE^ is self consistently generated because of plasma instabilities 
(e.g., ion-ion streaming instability (Papadopoulos et al., 1971) and is 
allowed to act directly on the ion motion, rather than through a model 
equation such as (l). Given this background, we now present the simulation 
results. 



The simulations were performed using a one dimensional ouasi-neutral 
hybrid code (Chodura, 1975; Winske and Leroy, 1984) whose description is 
given in Appendix I. The code incorporates ion kinetic effects by 
following the trajectories of "superparticles", numerically representing 
many actual ions, in 3 dimensions in velocity space and one dimension in 
configuration space. The electrons are approximated as a massless fluid 
described by momentum and energy transport equations; the electron pressure 
is assumed to be isotropic. These equations, with Maxwell's equations, are 
solved on a one dimensional spatial grid using moments calculated from the 
ion distributions. Anomalous resistivity, resulting in magnetic field 
diffusion and electron ohmic heating, is included as a parameter in the 
equations with a value consistent with that expected from current driven 
cross-field instabilities. 

The boundary conditions in the code correspond to the conditions we 
expect at early times. A dense ion beam is injected continuously into the 
simulation box which is initially filled with a tenuous stationary 
plasma. In our present results both ion populations are H + , but other ion 
species can be easily considered and will be presented elsewhere. We set 
up a finite width magnetic field compression within the leading edge of the 
ion beam. This represents the field compression expected from the cross- 
field beam motion which, because of finite resistivity, will diffuse into 
the beam. The width of the compression has been varied to study various 
physical phenomena. The magnetic field is held fixed at the edges of the 
simulation box. 

We present below the results of two simulations selected to illustrate 
the underlying physical coupling processes. The initial state is typical 
of the piston formation stage. It is shown in Fig. la, for the case of 
= V d /V a “ where V A is defined on the basis of the ambient (i.e., 
upstream) parameters. A magnetic pulse consistent with B/n = constant is 
introduced at the left hand boundary. In this paper we do not study the 
formation of the pulse but assume its nature in accordance with previous 
work mentioned in Section II, and its width L to Illustrate the important 
physics. In the figures the velocity is in terms of the upstream value of 
V A and of time in terms of the upstream value of Q *. The units of length 
are in terms of the ion cyclotron value computed for air with velocity 

5 






but with the value of the downstream (i.e., compressed) magnetic field. 
This representation has been selected for clarity of the underlying 
physics. We consider (1) a "broad" Lg = 2R^ (R^ in upstream values) pulse 
with M a 1 4 and n a /n a = 16 and (2) a "narrow" Lg = pulse with M A * 8 

and n a /n a * 64. Here V^, n d are the debris velocity and density, n a the 
air density and is defined on the basis of the upstream conditions. 

A. Broad Pulse (Lg * 2r Li ) 

This simulation is initialized with a magnetic pulse width Lg ~ 2r Li 
V d /V A = 4 and n d /n a = 16. We show a series of results for this run in Fig. 
1. Figure 1 displays four important quantities versus distance. The 
bottom curve represents the magnetic field magnitude, the middle curve (or 
"dots") represent particle velocities of the debris and air ions (the 

debris ions have V d ~ 4V^ while the air ions have V d ~ 0), and the top 
curve represents the electric field in the x direction. Figure la is at t 
= 0.0059 where is based on the upstream value of B. Note the broad 
pulse in B with sharp gradients at the leading and trailing edges. 
Associated with these gradients are strong laminar electric fields at these 
edges as shown in the top portion of Fig. la. These fields strongly affect 

the ions. Note in the middle panel that some debris ions are being 

accelerated ahead of the B pulse, and air ions are being picked up in a 

reflection process. 

Figure lb depicts the system at t * 0.0386 and several interesting 
features are developing. First, the B field pulse has broadened 

considerably. The leading edge of the pulse has been "dragged" out to 
x ~ 4r Ll , and the gradients are not as sharp. Second, the laminar electric 
field has fallen in intensity and is accompanied by a considerable amount 
of turbulence in the region 0 < r/r^ £ 2.4 caused by of the ion-ion 
streaming. This turbulence causes some momentum exchange between debris 
and air ions, and thermalization of the particles in this region. Finally, 
also at this time, a fraction of debris ions have been accelerated to a 
velocity greater than twice the initial velocity. 

Figure lc shows the system at t * 0.0874 The magnetic pulse 

continues to broaden and now extends to x ~ br^. A laminar electric field 
still exists at x ~ 3.6 r^. However, the turbulence in the region 0 < 







x/r^ < 3.6 has subsided since the relative velocity between debris and air 
ions is small. The ion motion in this region is dominated by the magnetic 
field and, in effect, is where "Larmor coupling" is taking place. Perhaps 
the most interesting region is 3.6 < x/r L1 < 6.0. It is clear that debris 
acceleration and air reflection are continuing. Furthermore, turbulence is 
developing at the debris-air interface (x ~ 4.0 r^). 

Finally, Fig. Id shows the system at t = 0.1542 The magnetic 

pulse is quite broad now, extending out to x ~ 8.0 r^. Very intense 
electrostatic turbulence is occurring at the debris-air mix in the region 
4.2 < x/r^^ < 8.2 which is producing coupling between the ion species and 
therraalizing the particles. Debris particles have been accelerated up to 

V ~ 3V d . In the region 2.0 < x/r^ < 4.0 the air ions have a velocity 

V ~ 5.0 V^ while the debris ions have a velocity V ~ 3.0 V^. In this 
region the momentum coupling is through the magnetic field (Larmor 
coupling) and the debris and air ions perform gyrations about the magnetic 
field. Very little turbulence is evident in the electric field, confirming 
the fact that v* < Sl a . 

Thus, in this simulation it is found that three important processes 
are affecting coupling, energization, and themaLlzatton: (1) laminar 
electric fields; (2) turbulent electric fields; and (3) magnetic field 
pick-up. The laminar fields are most prominent early in the run (Fig. la) 
and cause debris acceleration and air reflection. This is crucial to early 
time HANE processes since debris ions can be accelerated up to ~ 3 times 
their initial expansion velocity, and air ions are "picked-up". Turbulent 
ei.ecttic fields cause coupling and thermalization of the ion species. In 
this rjn they are most effective at the debris-air Interface, and ahead of 
the interface (Figs. 1c and Id). Finally, "Larmor coupling" is observed 
within a section of the broad magnetic pulse in Figs. Ic and Id, i.e., the 
debris and air ions have rotated about the B field so that the air ions 
have a larger velocity than that of the debris (V gx > V^ x ). This process 
supplements any incomplete part of the turbulent coupling and operates on a 


long time scal< 









B. Sari 


row pulse =0.1 r^) 

This simulation is initialized with a magnetic pulse width =0.1 
r^j and V^/V^ » 8 and n ( j/n a = 64. The results are shown in Fig. 2. In 

Fig. 2a we show the plasma system at t = 0.0206 . Already by this time 

a substantial number of air ions are being reflected from the laminar 
electric field at x ~ U5 r Li . Furthermore, debris ions have been 
accelerated up to V ~ 2.5 and extend out to x ~ 4r Li . In this region 
there is a considerable amount of electrostatic turbulence as seen in the 
electric field plot. Figure 2b is at t = 0.0411 0^*. The magnitude of the 
turbulence electric fields have increased from the previous time shown and 
are quite intense at the debris-air interface (x ~ 3.0 r Li ). The air tons 
are being "picked up" because of their turbulence and thermalized. The 
magnetic field has been dragged out to x ~ 7.0 r Li by accelerated debris 

ions. Finally, Figs. 2c and 2d show the ions and fields at t = 
0.0617 and t = 0.0771 respectively. The general features are 

similar to t = 0.0411 (Fig. 2b). There are intense turbulent fields at the 
debris-air interface which transfer momentum from the debris to air ions, 
and also thermalize the particles. Debris ions extend out to x > 10.0 r^ 
by t = 0.0617 0^, as well as the B field and low level turbulence. In 
this run there is no indication of Larmor coupling. One interesting 
feature seen in Fig. 2d is the "structuring" of the magnetic pulse. In 
particular, there is a sharp gradient that has developed at x ~ 4.8 r Li 
which produces a laminar electric field. This field acts to energize air 
ions to such that V ~ 1.8 at x ~ 5.0 r^. 

IV. DISCUSSION 

We have presented the results of simulations of very early time 
processes (i.e., t < 0^) using a ID hybrid code. The purpose of this work 
has been to elucidate important physical phenomena rather than to perform 
an actual early time simulation of a HANF. The focus of this study has 
been primarily on coupling mechanisms; an area that has been somewhat 
controversial within the DNA community. We have found that there is no 
"single" coupling mechanism but that several mechanisms can be effective at 
various times throughout the debris expansion phase. In the first 
simulation, initialized with a "broad" magnetic pulse, the initial coupling 
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occurs through laminar electric field at the leading edge of the pulse. 
Subsequent to this, turbulent coupling occurs at the initial debris-air 
interface while Larmor coupling occurs within the main body of the magnetic 
pulse. In the second simulation, initialized with a "narrow" magnetic 
pulse, the only coupling observed was because of the electric field 
(laminar and turbulent). 

A very important result of this study is the generation of very 
energetic debris ions. This is caused by the laminar electric field at the 
leading edge of the magnetic pulse. A fraction of the debris ions can be 
accelerated up to roughly three times the initial debris expansion velocity 
(i.e., ~ SVj). It is well-known that the initial formation of ion 
debris patches occurs on a time scale faster than that calculated from the 
initial expansion velocity. Thus, the accelerated debris ions observed in 
the simulations is consistent with HANE observations. 

Finally, this study dramatically highlights the importance of kinetic 
phenomena for the understanding of early time processes. Kinetic effects 
such as plasma turbulence, reflected and accelerated ions, and 
therraalization of the ions are not fluid processes, and cannot be recovered 
using only MHD codes. Detailed studies are necessary using codes such as 
the one described in this paper to properly understand the dynamics of 
early time HANE expansions. 
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Fig. 1 — Evolution of the ambient magnetic field (bottom curve), debris and air ions (middle 
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Fig. 1 (Corn’d) — Evolution of the ambient magnetic field (bottom curve), debris and air ions 
(middle curves/dots), and electric field in x-direction (top curve) as a function of time. The 
parameters are M A - VJ V A = 4, L a — 2 r Ll , and nj n a — 16. (a) t = 0.0059 ft, -1 ; (b) t = 0.0386 
nr 1 : (c) t = 0.0874 nr 1 : and (d) t = 0.1342 nr 1 . 
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nr 1 : (c) t =* 0.0874 nr 1 : and (d) t = 0.1342 nr 1 . 
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Fig. 1 (Corn’d) — Evolution of the ambient magnetic field (bottom curve), debris and air ions 
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parameters are M A >» VJ V A = 4, L B — 2 r Li , and njn a - 16. (a) t = 0.0059 0>) t * 0.0386 
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Fig. 2 (Corn’d) — Evolution of quantities described in Figure 1 but with Xf 4 =» VJ V 4 - 8, 
L B ~ 0.5 r u , and njn a - 64. (a) t = 0.0206 O' 1 ; (b) t = 0.0411 ftr 1 : (c) t - 0.0617 nr 1 : 
and (d> t - 0.0771 fl" 1 . 
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Fig. 2 (Cont'd) — Evolution of quantities described in Figure 1 but with \t 4 — VJ V A =- 8. 
L„ ~ 0.5 r Lh and n d /n a - 64. (a) t - 0.0206 flf 1 ; (b) t = 0.0411 Of 1 : (c) t - 0.0617 fi,"'; 
and (d) t - 0.0771 ft,' 1 - 
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APPENDIX 


HYBRID CODE DESCRIPTION 

We describe here the numerical details of the hybrid code used in this 
paper. This description is taken from Winske and Leroy (1983) with minor 
changes necessary due to the HANE application. Some of the techniques 
(e.g., particle-in-cell method) are standard and well known, while others 
(e.g., solution of the field equations) have been applied primarily to 
fusion-related problems and hence are generally not familiar to the space 
science community. The discussion is divided into five parts: (A) an 
overall description of the simulation methods, (B) treatment of the ions, 
(C) solutions of the field and electron fluid equations, (D) refinements to 
the basic model, and (E) discussion of the initial and boundary 
conditions. The methods pertaining to the first four parts are rather 
general and could thus be applied to any problem where a 'hybrid' (kinetic 
ion, fluid electron) description is appropriate. Specific adaptation of 
the method to the HANE study enters primarily through the last topic, the 
initial and boundary conditions. 

A. Overall Description 

The simulation method is particle-in-cell (Morse, 1970) with the 
electrons treated as a fluid (i.e., cell) quantity rather than as discrete 
particles. The assumption of a massless, charge neutral electron fluid 
thus eliminates the restriction of time and spatial lengths to he inverse 
electron plasma frequency and electron Debye length, usually associated 
with full particle codes. The hybrid models discussed here were originally 
developed for and extensively applied to dynamical studies of high density 
pinch experiments by Chodura (1975), Sgro and Nielson (1976) and Hamasaki 
et al. (1977). These models apply strictly to the case of only one spatial 
dimension; extension to two dimensions is possible, although nontrivial, as 
discussed later. Even with only one spatial dimension (which is along the 
x-axis and thus implies 3/3y = 3/3z * 0), all three velocity, magnetic 
field and electric field components are included. The simulation region 
has length L, divided into N cells, each of length Ax = L/N. There is an 
additional ('ghost') cell at each end of the system, which helps in keeping 
track of the particles entering and exiting the system and in setting up 
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the boundary conditions. The ceil quantities (i.e., fields, electron fluid 
properties, ion velocity moments) are specified at the cell centers. 
Because the restriction to electron spatial and temporal scales has been 
eliminated, much larger time steps can be used. Typically, the time 
step (At) is limited by the condition that ion gyromotion is well defined 
and the cell size is constrained by the condition that the fastest ions do 
not traverse one cell in one time step. The cell size and system length 
are also chosen to resolve length scales of interest to the problem. 

3. Dynamics 

The ion component is modeled by a discrete set of particles. The ion 
distribution is advanced in time by stepping forward each particle in time 
under the influence of the local, self-consistent Lorentz force. The 
motion of the ions in a four-dimensional phase space (v x , v y , v z , x) is 
solved by the particle-in-cell technique, using a second-order-accurate but 
non-reversible scheme (Nielson and Lewis, 1976). The equations for the ion 
advance (with superscripts denoting the time level) are: 

v° = v~ 1/2 + (h/2) E° 

v 1/2 - fv~ l/2 + h(E° + gB° + + P°) 


(Al) 


where h = At e/ra^, f * l-(h 2 /2)B^ • B^, g * fh/2)(v ^ 2 • B^), e is the 
charge and ra^ is the mass of the ions. The electric (e) and magnetic 
fields (b) are evaluated at the particle position x^. P is a mean friction 
force, P ■ -eg • J, exerted by the electrons as a macroscopic force 
only, J is the current and g represents a phenomenological anomalous 
resistivity which gives rise to Ohmic heating, as will be explained 
later. After all of the ions have been advanced, the ion density fn^) is 
updated, by averaging over the positions of all the particles in each 
cell. The ion mean velocities ) are also needed at time level 1. The 
velocities are pushed ahead one half time step, 

v l/2 x B l 

V 1 =■ V 1/2 + Ch/2)fE l + :: --—— + p 1 ), (A2) 
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using the fields at the new particle positions x\ in order to accumulate 
the velocity moments V^. 

The assumption of charge neutrality then implies 

n = n. = n, (A3) 

e i 


and the continuity equation gives 


at 


K + »,) ' 


3 _ 

3x 


en fv 

ix 


V 1 
ex' 


(AA) 


V 


ex 


(A5) 


Thus, the electron density and one component of the electron fluid velocity 
are determined. 


C. Field-Fluid Equations 

With the assumption of zero electron mass the electron momentum 
equation reduces to 


dV 

~e 

n m -j— 
e e dt 


) “ 3x J 


(A6) 


where p e is the scalar pressure and is the velocity of the electron 
fluid. The transverse (y and z) components of (A6) can be written as Ohm's 
law 


J 



(A7) 


In general, the magnetic field is inclined at some angle with respect to 

the x axis: 3 = |b|( b , b , b ) with b^ + b^ + b^ = 1, and if a-fo.) 

~ 1 1 ’ x’ y’ z' x y z Hi- 

represents the conductivity parallel (perpendicular) to the magnetic field, 
the conductivity tensor is (Krall and Trivelplece, 1973): 
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-i 1 

2 = 3 " ° y 

a 

a a 

y z 

a , 

2 x 

a x a 3 


(A8) 

• ~ 

: 





• .. 


where 





i 

*1 = a i + 

K - 'A 1 , 



• 


a 2 = °1 + 

- ° a >y b y 





a 3 = *L + 

t". - °i) b z b z 



• . 


% • K- 

°lVz 




b 

a = fa. - 
y v 1 

°A b z 



• . 


a Z = K " 

a, ]b b . 

1' x y 





The transverse fields 

can be expressed 

in terms of 

the vector 



potentials in the usual manner f§, “ 2 x A; E * - 

fl/c)f3A/3t)1: 


¥~*~ 


3A 

\ " JT 






3A 

B z " 3^ 

E - - I 

y c 3t 



(A9) 



. 3A 

E - - I —2 

z c 3t 




• 


and 2*3=0 implies 






= constant. 



(A10) 

* 


(The last field component, E 

x , is derived from 

the x-component 

of (A6) as 



will be shown later.) 
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Neglecting the displacement current. Ampere's L 


can be expressed 



Substituting in Ohm's law (A7) and the definitions of the 
components, Eq. (11) becomes a coupled system of equations: 

3“A v 3A v 3A z 3A 3A z 

TT = f l ’ lt~ ’ 3x ’ 

3x 

3 2 A 3A 3A 3A 3A 

_5. = f _5. _5.) 

9 2 r 2 '3t ’ 3t ’ 3x ’ 3x J 


In order to solve this system a fully implicit, space ceni 
difference scheme is used (Richtmyer and Morton, 1967). If we let u^ 
at the jth cell position and n-th time level (and w" = A z in the same : 
the left hand side of (A12) is differenced as follows: 


n+l 

Vl 


u n+l 


3x 


(A13) 


Ax 


while on the right hand side 

n+l n 


3A 


3t 


At 


n+l 

Vl 


n+l 

U i-1 


When this is substituted into (A12), the resulting system of equations can 
be written symbolically as 


* X44.1 


X, + C. 


X. 


(, 






X 


•j 


n+1 

U j 

n+1 


(A16) 


and Aj, B^, and Cj are 2x2 matrices that depend only on components of a and 
(evaluated at cell j at time level n) and constants, while D also 
depends on u?, w^. Thus, all the coefficients A, B, C and D are explicitly 
known. Note that there are N cells in the computation mesh j*l, 2, ... N 

with ghost cells (.1*0 and 1=N+1) at each end. 


The set of equations (A15) is easily solved. Assuming a solution of 
the form 


Sj ■ tj • Vi + T S j - 0 , 1, ... s am 

and substituting it into (A15), it then follows that 

-- i 8j -ij-.r 1 8j 

£j - !lj ♦ Sj • V,r l • fBj - Sj • <* 18 > 

Because all the A., B.., and are known, if E^ and Fq are known (from 
boundary conditions, as explained Later), (A18) can be solved in ascending 
order to obtain E^, F ^; Eg, Fg, •••, g N , F^. Then if X^ can be 
determined (again from boundary conditions), (A17) can be solved in 
descending order to obtain X^, X^j, • ••, Xq* 

Once the components of the vector potential are known, the new values 
of the electric and magnetic field follow from (A9) and the components of 
the current from (All), again using the differencing in (A13) and (A14). 
The transverse components of V are then easily found from the current and 
the ion velocity moments: 

V e.j = V ij " J j /en .1 = y or z. (A19) 
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The last electron fluid quantity to be calculated is the temperature. 
The electron energy equation can be written as 


2 _ 3 _ 
2 3t 


n T V 
e e ex 


3V 


ex 

3x 


Q. 


(A20) 


2 

The source term Q includes resistive heating (nJ ) and loss mechanisms, 
such as thermal conduction or radiation, depending on the application. 
Letting 


(A20) can be rewritten in the form (with y * 5/3) 
. , 3V 


(A22) 


The differencing for 3p g /3t and 3V ex /3x is identical to (A14); for 
stability reasons "donor cell" differencing is used for the convective terra 
(Richtmyer and Morton, 1967): 


3p V . 
e m exj 
c 3x Ax 


n+l 

’ej 

n+l > 

' P ej-2 ^ 

V exj 

n+l 

e.1+1 

n+l 

• Vi > 

V exj 


again, a tridiagonal system is obtained, 
A jXj+ i + B jXj + 


> 0 

< 0 


(A23) 


(A24) 


except that now ail quantities are simple scalars. As before, an assumed 
solution of the form 


Ej - - [Bj + c jVl r A, 


(A25) 


24 






j 


(A26) 


* c»j + cj'j-i r 1 (»j - CjPj,, ]. 

Appropriate boundary conditions lead to a determination of Eg, Fg (and then 
to all Ej, Fj in ascending order, using (A26)) and X N+1 ( and thus to all Xj 
in descending order through (A25)). Since n e (= n^ is known, T g is then 
obtained from (A21). An alternative to (A22) which is sometimes useful 
(Sgro, 1978) is to write the energy equation (A20) in terms of the 
entropy [s = in(p e n e ^)l instead of the pressure. 

Finally, the x component of the electron momentum equation (A6) can be 
solved for E x : 


■ Cv 


ne 3x 


(A27) 


This electric field is needed to maintain charge neutrality. 

The computation loop through one time step can thus be summarized as 
follows: 

1. Advance the ions one time step and then calculate the ion velocity 
moments (n^ ■ n fi , V ix 3 V ex» v iy* V iz ) * 

2. Compute the plasma conductivity g, according to some prescription. 
(The resistivity g is g ). 

3. Solve the coupled equations for A y , A z ; we can then easily 
calculate E y , E z , B y , B z (B x ■ constant). From the currents 
(Ampere's law) and the ion moments calculate the other two 
electron velocity components, V gy and V ez . 

4. The electron temperature is obtained next from the solution of the 
differenced energy equation. 

5. E x is then calculated from the x-component of the electron 
momentum eauation. 

6. Thus, all field and electron fluid components are known and we are 
ready to move the ions again. 
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D. Refinements 


The simulation scheme described thus far Is rather general, the only 
assumptions being quasineutrality, zero electron mass and one spatial 
dimension. All of these conditions can be replaced by more appropriate 
ones, as the physical situation dictates. For example, for the study of 
low frequency ion waves an adiabatic electron model is more appropriate; 
such a model has been successfully used by Okuda et al. (1978). In some 
situations phenomena in the lower hybrid frequency range are of interest. 
In this case electron inertia effects are non-negligible; they have been 
included in the one-dimensional model of Liewer (A1976) and the two- 
dimensional model of Hewett and Nielson (A1978). Two-dimensional 
simulation models with m g = 0 have been successfully employed by Byers et 
al. (1978), Hewett (1980) and Harned (1982). 

The HANE simulations described in this paper take the conductivity to 
be constant. A further level of sophistication is to include anomalous 
processes due to the microphysics, which is occurring on time and distance 
scales shorter than those resolved in the hybrid model, by means of more 
complicated transport coefficients. This is done by expressing the 
conductivity (or resistivity) as a sura of two terms, one representing 
classical (Spitzer) effects and the other due to anomalous effects arising 
from microinstabilities due to cross-field currents. Two types of 
anomalous transport coefficients have been used successfully to model the 
behavior of laboratory plasmas: one type uses a serai-empirical expression 
(Chodura, 1975; Sgro and Nielson, 1976), the other is based on a 
quasillnear analysis of known instabilities (Davidson and Krall, 1977). 
The transport processes can be further refined to include electron thermal 
conductivity (Sgro, 1980) and anomalous ion heating (Hamasaki et al., 
1977). Multispecies ions can also be included in a straightforward manner 
(Sgro, 1980; Sgro and Winske, 1981). 
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E. Initial and Boundary Conditions 

In order to simulate the piston formation stage of the HANE problem, 
appropriate initial and boundary conditions have been implemented in the 
hybrid code. We assume initially that, a dense cool debris plasma flows 
across the magnetic field into a tenuous stationary air plasma. In this 
paper, we assume the debris and air to consist of the same ion species. 
This assumption is not necessary and we are presently investigating the 
effects of different debris composition. We assume further that there is a 
magnetic compression of finite width imbedded in the leading edge of the 
debris stream, formed as a result of the preceding magnetic compression 
stage of the coupling is simply maintained at the right boundary. Ions of 
all species are effectively absorbed at the boundaries. The final boundary 
condition is on the magnetic field which is held constant at the ambient 
(air) value at both boundaries. The magnetic compression is taken to be 
proportional to the electron compression. The width of the magnetic 
compression is an input parameter in the code. Thus initially in the 
simulation space we establish a debris ion stream extending a finite 
distance from the left boundary into the uniformly distributed air ions. A 
step function compression in the magnetic field, narrower than the debris 
stream, is initialized at the leading edge of the debris. The magnetic 
field is initially uniform on both sides of the compression. During the 
simulation runs, debris ions are continuously injected from the left 
boundary. This is accomplished by re-initializing the debris ion 
distribution in the left ghost cell each time step. 
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WASHINGTON, D.C. 20545 
01CY ATTN DOC CON DR 




50VERNMENT 


JTE FOR TELECOM SCIENCES 

•iAL TELECOMMNICATIONS & INFO 

IN 

ER, CO 80303 

t ATTN A. JEAN CUNCLASS ONLY) 
1 ATTN W. UTLAUT 
1 ATTN D. CROMBIE 
1 ATTN L. BERRY 


CHARLES STARK DRAPER LAB 
S 55 TECHNOLOGY SQUARE 

CAMBRIDGE, MA 02139 
01CY ATTN O.B. COX 
01 C Y ATTN J.P. GILMOP 


COMSAT LABORATORIES 
LINTHICUM ROAD 
CLARKSBURG, MD 20734 
01C Y ATTN G. HYDE 


YAL OCEANIC & ATMOSPHERIC ADMIN 
ONMENTAL RESEARCH LABORATORIES 
TMENT OF COMMERCE 

ER, CO 80302 

Y ATTN R. GRUBB 

Y ATTN AERONOMY LAS G. REID 


:TMENT OF DEFENSE CONTRACTORS 


CORNELL UNIVERSITY 
DEPARTMENT OF ELECTRICAl 
ITHACA, NY 14850 

01CY ATTN D.T. FARLE' 

ELECTROSPACE SYSTEMS, II 
BOX 1359 

RICHARDSON, TX 75080 
01CY ATTN H. LOGSTON 
01C Y ATTN SECURITY (F 


IPACE CORPORATION 
BOX 92957 
kNGELES, CA 90009 
:Y ATTN I. GARFUNKEL 
lY ATTN T. SALMI 
:Y ATTN V. JOSEPHSON 
:Y ATTN S. BOWER 
CY ATTN D. OLSEN 

CTICAL SYSTEMS ENGINEERING CORP 
> CONCORD ROAD 
INGTON, MA 01803 
CY ATTN RADIO SCIENCES 

IN RESEARCH ASSOC., INC. 

RUTLAND DRIVE 
IN, TX 78758 
CY ATTN L. SLOAN 
CY ATTN R. THOMPSON 


EOS TECHNOLOGIES, INC. 
606 Uilshire Blvd. 

Santa Monica, Calif 90i 
0ICY ATTN C.B. GABBAf 
01CY ATTN R. LELEVIEf 

ESL, INC. 

495 JAVA DRIVE 
SUNNYVALE, CA 94086 
01CY ATTN J. ROBERTS 
01CY ATTN JAMES MARSi 

GENERAL ELECTRIC C0MPAN1 
SPACE DIVISION 
VALLEY FORGE SPACE CENTI 
GODDARD BLVD KING OF PRl 
P.O. BOX 8555 
PHILADELPHIA, PA 19101 
01CY ATTN M.H. BORTNI 
SPACE SCI LAB 


iLEY RESEARCH ASSOCIATES, INC. 
BOX 983 

: I f Y rft 94701 


GENERAL ELECTRIC COMPAN 









107 COBLE HALL 

150 DAVENPORT HOUSE 

CHAMPAIGN, IL 61820 

(ALL CORRES ATTN DAN MCCLELLAND) 
01CY ATTN K. YEH 

INSTITUTE FOR DEFENSE ANALYSES 
1801 NO. BEAUREGARD STREET 
ALEXANDRIA, VA 22311 
01CY ATTN J.M. AEIN 
01C Y ATTN ERNEST BAUER 
01CY ATTN HANS WOLFARD 
01CY ATTN JOEL BENGSTON 

INTL TEL S TELEGRAPH CORPORATION 
500 WASHINGTON AVENUE 
NUTLEY, NJ 07110 

01CY ATTN TECHNICAL LIBRARY 

JAYCOR 

11011 TORREYANA ROAD 

P.O. SOX 85154 

SAN DIEGO, CA 92138 

01CY ATTN J.L. SPERLING 


P.O. BOX 504 
SUNNYVALE, CA 94088 
01CY ATTN DEPT 60-12 
01CY ATTN D.R. CHURCHILL 

LOCKHEED MISSILES 8 SPACE CO., INC. 
3251 HANOVER STREET 
PALO ALTO, CA 94304 

01CY ATTN MARTIN WALT DEPT 52-12 
01CY ATTN W.L. IMHOF DEPT 52-12 
01CY ATTN RICHARD G. JOHNSON 
DEPT 52-12 

01CY ATTN J.B. CLAD IS DEPT 52-12 

MARTIN MARIETTA CORP 
ORLANDO DIVISION 
~.U. BOX 5837 
ORLANDO, FL 32805 

01CY ATTN R. HEFFNER 

M.I.T. LINCOLN LABORATORY 
P.O. BOX 73 
LEXINGTON, KA 02173 

01 CY ATTN DAVID M. TOUI F 






735 STATE STREET 
SANTA BARBARA, CA 93101 
01CY ATTN P. FISCHER 
01CY ATTN U.F. CREVIER 
01CY ATTN STEVEN L. GUTSCHE 
01C Y ATTN R. BOGUSCH 
01CY ATTN R. HENDRICK 
01 C Y ATTN RALPH KILB 
01 CY ATTN DAVE SOWLE 
01C Y ATTN F. FAJEN 
01CY ATTN M. SCHEIBE 
01CY ATTN CONRAD L. LONGMIRE 
01CY ATTN B. WHITE 
01CY ATTN R. STAGAT 

MISSION RESEARCH CORP. 

1720 RANDOLPH ROAD, S.E. 
ALBUQUERQUE, NEW MEXICO 87106 
01C Y R. STELLINGWERF 
01CY M. ALME 
01CY L. WRIGHT 

MITRE CORPORATION, THE 
P.O. BOX 208 
BEDFORD, MA 01730 

01 CY ATTN JOHN MORGANSTERN 
0 1CY ATTN G. HARDING 
01C Y ATTN C.E. CALLAHAN 

MITRE CORP 

WESTGATE RESEARCH PARK 
1820 DOLLY MADISON BLVD 
MCLEAN, VA 22101 
01CY ATTN W. HALL 
01CY ATTN W. FOSTER 

PACIFIC-SIERRA RESEARCH CORP 
12 3 AO SANTA MONICA BLVD. 

LOS ANGELES. CA 90025 


P.O. BOX 3027 
BELLEVUE, WA 98009 

01CY ATTN E.J. FREMO 

PHYSICAL DYNAMICS, INC. 
P.O. BOX 10367 
OAKLAND, CA 94610 
ATTN A. THOMSON 

R ft D ASSOCIATES 
P.O. 90X 9695 
MARINA DEL REY, CA 902 
01CY ATTN FORREST GI 
01CY ATTN WILLIAM B. 
01CY ATTN WILLIAM J. 
01 CY ATTN H. ORY 
01CY ATTN C. MACOONA 
01CY ATTN R. TURCO 
01CY ATTN L. DeRAND 
01CY ATTN W. TSAI 

RAND CORPORATION, THE 
1700 MAIN STREET 
SANTA MONICA, CA 90406 
01 CY ATTN CULLEN CRA 
01CY ATTN ED BEDROZI 

RAYTHEON CO. 

528 BOSTON POST ROAD 
SUDBURY, MA 01776 

01 CY ATTN BARBARA AO 

RIVERSIDE RESEARCH INST 
330 WEST 42nd STREET 
NEW YORK, NY 10036 

01CY ATTN VINCE TRAF 














END 

FILMED 

7-85 


i 

DTIC i 



